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Abstract 

O 

The intranuclear cascade model overestimates the multiplicity of pro- 
| duced mesons in nucleus-nucleus interactions without taking into ac- 

^ ' count meson and baryon resonance production. Inclusion of the reso- 

nances leads to decreasing multiplicity of mesons, neutrons and protons. 
■ In order to overcome the problem, it is proposed to use the FRITIOF 

model adapted to low energies in a combination with the reggeon theory 
\ inspired model of nuclear destruction. It is shown that the combina- 

' tion allows one to reproduce satisfactory the meson and baryon yields 

in the nucleus-nucleus collisions at the energy of 3.3 GeV/nucleon. The 
£^ . combined model works faster then typical quantum molecular dynamic 

model, and allows one to estimate the data needed for creation of electro- 
g ; nuclear amplifier. 
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Data on nuclear reactions at the energies of hundred MeV and GeV are required for 
multiple purposes such as long-lived radioactive waste transmutation, material analysis, 
nuclear medicine as well as research of the cosmic ray effects on spaceships and astronauts. 
Experiments on measuring the data are costly to be carried out and there is a limited 
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number of facilities to make them. Therefore, reliable computer models for a simulation 



of the reactions are created to provide the necessary data. Most of them are using the 
ideas of the cascade-evaporation model (CEM) (see 01 - ]§])• 

The cascading of nucleons and 7r-mesons was considered only by the first variants of 
the model. A good description of hadron-nucleus interactions was reached at the studies. 
The best result was obtained by S. Mashnik 0. Though, the application of the model 
for description of nucleus-nucleus collisions has shown that the model gives a satisfactory 
yield of the nucleon but overestimates the meson production. Taking into account meson 
and baryon resonances production is one of the possible ways to solve the problem. 

A number of authors were trying to do this, and they usually obtained a decreas- 
ing yield of mesons and baryons. It is natural because the effective decreasing of the 
multiplicity of the produced particles leads to a less powerful cascading. So, a problem 
of a self consistent description of the meson and baryon yields in hadron-nucleus and 
nucleus-nucleus interactions appeared. 

Let us note that the nuclear destruction mechanism and the procedure of the excitation 
energy calculation were not changed in the mentioned approaches. Maybe, they led to 
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the unsatisfactory results. In this paper we consider a synthesis of the FRITIOF model 
|| [J which takes into account the resonances production and the reggeon theory inspired 



model of the nuclear destruction 10, 11 



The FRITIOF code |J, that is a program of Monte Carlo simulation of the inelastic 
hadron-hadron, hadron-nucleus and nucleus-nucleus interactions, is very popular in high 
energy experimental physics. It is explained by its access, its physical ideas simplicity 
and a defined beauty of the code by itself. It is easy to use. 

The FRITIOF model || || assumes that an excitation of hadrons into continuum mass 
spectra takes place in the inelastic hadron-hadron collisions. In the case of the hadron- 
nucleus or the nucleus-nucleus interactions the excited hadrons can suffer an additional 
collisions with the nuclear nucleons and go into more excited states, or de-excite. The 
excited hadrons are considered as quark strings, and the corresponding quark model 



12, 13] is used for description of their decay. The probabilities of the multiple collisions 



are calculated within the Glauber model |L4| |19[. The inelastic collisions are usually 



considered only [|T^]. In order to reproduce the baryon yield, we introduce elastic re- 
scattering too. 

It is assumed that the program can not be used at the relatively low energies as the 
hypothesis about the creation and decay of the quark strings is not valid. Attempts 
of compelling the program to operate at the energies below 5-10 GeV/nucleon for AA- 
interactions usually failed. Though the analysis of the code operation shows that the 
program cycles due to its quite simple-hearted interpretation of the Fermi-motion of the 
nucleons. A change of the Fermi-motion simulation algorithm [^] allowed one to decrease 
a formal limit of the model application region. 

Now in order to use the model at the intermediate energies, one needs to correct a 
scheme of the fragmentation of the quark string with low masses because in this case the 
excited hadrons have low masses, too. Below the corresponding changes of the model will 
be presented which allow one to reach a defined success. To solve the problem, we have 
used the experimental data on the neutron-proton interactions at momentum 1.25 - 5.1 
GeV/c given amiable by the Yu. A. Troyan's group of the neutron-proton interaction 
study. 

One of the disadvantages of the FRITIOF code is an omission of the slow particle 
cascading into the nuclei. Under the "cascading" one usually understands the standard 
intranuclear cascade scenario (see, e.g., (l[|), which neglects the quantum mechanical ef- 
fects. The quantum mechanical description of the particles cascading in the nuclei can be 
achieved in the framework of the reggeon theory. According to the reggeon approach |22 



a consideration of the cascade interactions necessitates a calculation of the yields of the 
so-called enhanced diagrams to the elastic scattering amplitude. Using the procedure for 
the calculation proposed in [[23| and Abramovski-Gribov-Kancheli cutting rules one 
can obtain the positive defined cross-sections of the inelastic processes. At first glance, 
this returns one to the classical cascade picture of the interactions. However, there is 
an essential difference. The cascade model assumes that the cascade is developed in a 
three-dimensional space of a target nucleus. According to the reggeon approach |I| , the 
"cascade" of the reggeon exchanges occurs in a two-dimensional space of projected radius- 
vectors of nucleons on the plane perpendicular to the momentum of projectile particle (on 
a plane of impact parameter), with a "cascade power" independent of the multiplicity of 
the produced particles and defined by the reggeon vertex constants and the size of nu- 
cleus. We shell give a corresponding algorithm in Sec. 3. The main calculation results 
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are presented in Sec. 4. We are starting with a short description of the main assumptions 
of the FRITIOF model. 



2 Theses of the FRITIOF model 

The FRITIOF model assumes that the hadron-hadron interactions are of two-particle 
character 

a + b — >a' + b', (1) 

where a' and b' are a and b hadrons in excited states. The kinematics of the reaction is 
determined as follows: in the center of mass of colliding hadrons the energy-momentum 
conservation low has the form: 

E a + E b = E a > + Ey = 

Paz+Pbz = Pa'z+Pb'z = 0, (2) 
= Pa<±+Pb>±, 

where E a and E b (E a >,E b /)- energies of the initial (final) hadron a and b (a',b r ), p az and 
Pbz - longitudinal momentum components (projected momenta on the interaction axis). 
Adding and subtracting the first two equations from (Q), we get 

P+ + P+ = P+ + P+ 

P- + P b - = P-,+Py (3) 
= Pa'±+Pb'±, 

where P + = E + p z , P~ = E — p z . 

At high energies: P^ ~ m a'/^ I Pa'z \, Py — m 6'/2 I Pb'z \ ■ Thus, the P~, and P b t 
distributions used in the code have the form 

dW ~ dP-,/p-,~dml,/ml„ 

dW ~ dP b t/P b t ~ dml/ml. (4) 

The limits of P~, and P b t are defined as 

[P-,P b -], [P b + ,P+]. (5) 

The distributions (Q) are typical for the so-called high-mass diffraction dissociation pro- 
cesses. 

Knowing P~, , P b t, p a >±, Pb'X and determining Pjr, P b 7 from Eqs. (|| ), one can find the 
masses of the excited hadrons a' and b' . p a f± and py±_ are sampled according to the low 

dW ~ exp{-f a ,J < p 2 ± >) d 2 p a r X . (6) 

In case of hadron-nucleons interactions the kinematics governed by Eqs. @, @ 
is applied to the first collision of the projectile nucleon with one of the target nucleons 
(a + N% — > a' + N[). For the second collision (a' + N 2 — > a" + A^), analogous relations are 
used, but (||) is replaced for 

[P 7,P^J, [P^ 2 ,P Q t]. (7) 
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As a result, the consequent collisions involve a systematic increasing of the mass of hadron 
a if transfers of the transverse momentum are small. 

A similar approach is also applied to simulate the nucleus-nucleus interactions. Here 
the reactions a' + V —>■ a" + b" are considered. The above distributions on P~, and P 6 t are 
replaced by those on P~„ and P b t and the limits of P~„ and P b t are redefined as 

[P-,,Py], [P 6 t,P t]. (8) 

At relatively lower energies of the order of 5-10 GeV/nucleon, the FRITIOF model 
without taking into account hadron de-excitation overestimates the multiplicity of the 
produced particles in hA- and AA-interactions. To this end we have changed the condition 
(|8|) by the following giving an allowance for considering the excitation process with the 
increasing mass and the de-excitation process with decreasing mass: 

[Pa, \f^>v - m b], [P& + , \/~^>v ~ m a \. (9) 

P~ and P b + are calculated at s a b = s a >b>, m a i = m a , my = m^. The minimal values of P~ 
and P b + , being the lower limits in ([5|), are obviously reached in the reaction a' + b' — > a + b. 
The maximal values are achieved at p a » z = py z — when hadrons come to rest in the c. 
m. frame of iViV-collision. 



The reactions (H), or 



a + b — >a' + b, (10) 
a + b — >a + b', (11) 



are the so-called diffraction dissociation reactions. The reactions (|T0|), ( |TTD are one- vertex 
diffractions, the reaction ([!]) is a double vertex diffraction. It is obvious that a minimal 
mass of the excited nucleon in NN-interactions can be equal to m a i = m^ + = 1080 
MeV. In the FRITIOF model the minimal mass is equal to 1.2 GeV. Due to this, the 
one vertex diffraction with excitation only one hadron is possible at ^snn < 2.4 GeV, 
or at Piat, < 1.91 GeV/c in NN-interactions. There can be a two vertex diffraction at 
higher energies. A relation between the cross-sections of the processes is determined by 
the hadrons a' and b' mass distributions. 

The excited hadrons a! and b' are considered as quark strings, and the corresponding 



quark model is used for the simulation of their decay |12|, [13| . It is assumed that the 
quark model can be used at sufficient large string masses what can be created at high 
energies. Thus the FRITIOF model was used mainly at high energies. One can expect 
that the FRITIOF model predictions fall into a contradiction with experimental data 
on the processes where the states with low masses can appear. In order to study the 
situation, we turn to the data on 7r~-meson and proton distributions in rip-interactions 
at P n = 1.25-5.1 GeV/c 

Fig. |I] shows the experimental and calculated 7r~-meson rapidity distributions, y = 
| \n(E + P z )/(E — P z ), where E and P z are laboratory energy and longitudinal momentum 
of 7r~-meson, respectively. As seen, the distributions calculated according to the original 
code (dashed lines) are of two-bump structure more pronounced at low energies. It seems 
like it is a circumstance of the assumed diffractive character of the interactions: the bump 
at large rapidities is caused by a projectile hadron diffraction, the bump at low rapidities 
is connected with a target nucleon diffraction. Though, at the neutron momentum of 
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1.25 GeV/c a diffraction system with mass of 1.2 GeV (the minimal mass of excited 
nucleon assumed by the FRITIOF model) must be in a rest in the center of mass system, 
and there must not be a subdivision of the fragmentation regions. Thus, we conclude 
that the two-bump structure of the calculated distributions is not a circumstance of the 
diffraction dissociation of the hadrons. It only reflects an anysotropy of the low mass 
string decay. A direct simulation of the decay of the strings with low masses has shown 
that the bi-module structure is typical for the decay of the strings with masses lower 
than 1.7 GeV. Two-particle channel is dominating in the decay of such strings. At higher 
masses the multi-particle channel gets more probable, and the calculated distributions get 
more regular. 

Taking into account the character of the experimental distributions at P n = 1.25, 1.73 
GeV/c, it seems reasonable to simulate an isotropic decay of the strings with low masses 
in the case of the two-particle decay channel. A boundary value of the string mass of 1.7 
GeV below of which we have used the proposed procedure, was chosen requiring a good 
description of the data. As seen in the Fig. |I| (solid curves), this allows us to describe the 
experimental data quite well. 

A more complete situation appears in the description of the 7r~-meson distributions 
in transverse momentum, P?- The original model predicted the average transverse mo- 
mentum larger than the experimental ones (see the dashed curves in Fig. |2|), though we 
had changed the character of the low mass string decay. Since at P„ = 1.25, we simu- 
late an isotropic decay of the excited hadrons, the maximum of the distribution in Pt 
is determined by mass of the decayed system. Thus, in order to decrease the average 
transverse momentum, the minimal mass of the excited nucleon state was decreased to 
the value of 1.1 GeV. This, as seen from Fig. ||, gave a better result at P n < 2 GeV/c. 
At larger energies one needs to take into account the transverse momentum transferred 
by the colliding nucleons. 

An analysis of proton spectra gives some additional information, and allows one to 
determinate the model parameters more exactly. For the analysis we have used the data 
on the following reactions 

np — > pp7i~; (12) 
np — > ppir~n ; (13) 
np — > np7T + 7T~. (14) 

In the reaction (|12"D according to the model the diffraction dissociation processes of 
the projectile particle is a dominant one. In the reaction ( |T3"D the two-vertex diffraction 
processes do the same. At last, in the reaction ( |Hj) we have equal yields of the one- 
vertex and two-vertex diffraction. Thus, the study of the reactions allows one to check 
the different components of the model. 

Fig. |3| gives the experimental and calculated rapidity distributions of the protons. 
The experimental data are presented by the histograms, the original model calculations 
- by the dashed curves, and the last calculations - by the solid curves. As seen, there is 
a two-bump structure of the experimental distributions on the reactions (0). The bump 
at small rapidities is caused by the saved target protons. The bump at large rapidities is 
connected with the protons created in the projectile particle diffraction. The distribution 
of the reaction (|T3|) has no structure. At last, in the reactions (|T4T ) there is a dominant 
production of the protons in the target fragmentation region. 
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The original model calculations are in agreement with the data on the reaction (|I4"D. 
At the same time, the calculated distribution for the reaction (|i~2"D has a bump in the 
region y ~ 1.7 what does not observe at the experiment. For its elimination we introduce 
a charge exchange between the nucleons in 50 % of the two- vertex diffraction. This allowed 
us to improve the proton spectra description in part. 

Enumeration of the changes made in the FRITIOF code 



1. The minimal mass of the excited nucleon decreases from 1.2 GeV to 1.1 GeV; 

2. In case of the two-particle decay channel of a string with a mass lower than 1.7 GeV 
the isotropic decay is simulated in the center of the mass system; 

3. The charge exchange between the colliding nucleons is allowed in 50 % of the two- 
vertex diffraction; 

4. The value of the average square of the transverse momentum which transferred 
between the colliding nucleons increases from 0.08 (GeV/c) 2 to 0.15 (GeV/c) 2 . 

3 Simulation of nuclear destruction at the fast stage 
of interaction. 



3.1 Determination of the number of knocked-out nucleons 



In the last few years there have been some successful attempts to describe the hadron- 
hadron elastic scattering at low and intermediate energies (below 1 — 2 GeV) within the 
quark-gluon approach (see Refs. [^5] 
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In Ref. [^5| the amplitudes of nn— , 
Kit— and NN— scattering were found and an agreement of the theoretical calculations 
with corresponding experimental data was reached at the assumption that in the elas- 
tic hadron scattering one-gluon exchange with the following quark interchange between 
hadrons takes place (see Fig. |]a). At high energies two-gluon exchange appropriation 
(Fig.[|b) works quite well (see Ref. [p9] , p0| and [[H]]). What kind of exchanges can 



dominate in hadron-nucleus and nucleus-nucleus interactions? 

The simplest possible diagrams of the processes with three nucleons are given in Fig. 
[5]. Calculation of their amplitudes according to Refs. [p5]-[p8 is a serious mathematical 
problem. It can be simplified if one takes into account an analogy between the quark- 
gluon diagrams and the reggeon diagrams: the quark diagram of Fig. |]a corresponds 
to a one-nonvacuum-reggeon exchange diagram; the diagram of Fig. |]b describes the 
pomeron exchange in the t— channel; the diagram of Fig. [|a is in a correspondence with 
the enhanced reggeon diagram of the pomeron splitting into two non-vacuum reggeons. 
The three pomeron diagram (Fig. ||d) represents a more complicated process. It is rather 
hard to find a correspondence between the reggeon diagrams and the diagrams of Fig. 
|b,|c. 

The reggeon parameters and the functional forms of the amplitudes of 3-reggeon pro- 
cesses are well known. The constants of the reggeon interaction vertexes are poor deter- 
mined. The 3-pomeron vertex constant Gppp is well established (Gppp = 1.35~ 2 (GeV^) 2 , 
Ref. |23]|). There are only old data [32| and the estimations of Ref. on the values 
of other constants - Gppp and Gppp, which are large. Nevertheless, we believe that the 
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properties of the reggeon amplitudes must be taken into account in consideration of the 
nuclear destruction. 

It is obvious that the processes like that in Fig. [5]d can not dominate in the elastic 
hadron-nucleus scattering because they are accompanied by production of a high mass 
diffraction beam of the particles in the intermediate state. Thus, their yields are dumped 
by a nuclear form-factor. According to the same reason, the yields of the processes like 
ones in Figs. |5|a, [5|b can be small, too. If it is not so, one will expect large corrections to 
Glauber's cross-sections. The practice shows that the corrections to the hadron-nucleus 
cross-sections must be lower than 5 - 7 %. 

The yield of the diagram of Fig. |5]c gives a correction to Glauber's one-scattering 
amplitude. There must be analogous corrections to the other terms of Glauber's series. 
The sum of the corrections must lead to small effects in the elastic small angle scattering 
because the corrections are large at small impact parameters. So, they can manifest 
themselves at large scattering angles. We assume that they have a big influence on the 
inelastic process characteristics, too. 

According to the reggeon theory, a description of the inelastic reactions can be reached 
in a consideration of the different cuts of the reggeon diagrams. Here the Abramovski 
- Gribov - Kancheli cutting rules are often used. The corrections to them were 



discussed in Ref. [22] in the application to the problem of a particle cascading on the 



nucleus. As was shown in Ref. [p2 |, summation of the yields of enhanced diagrams allows 



one to describe increasing of the one-particle spectra in the target fragmentation region. 
At the same time, the authors of Ref. p2| did not take into account the shadowing effects 



considered in Ref. pi ] 



Here we have to note that the yields of the diagrams like that shown in Fig. |5|c have 
no shadowing corrections. The yield of the enhanced diagram of Fig. |5]a has a form 

Y a ~ exp[-(h - b 2 ) 2 /3r 2 a - (h - b 3 ) 2 /3r 2 a - (b 2 - b 3 ) 2 /3r 2 a ] 

where bi,b 2 and 63 are the impact coordinates of the nucleons. At the same time, the 
yield of the diagram of Fig. [5]c according to Refs. [25] - |28| is given by 

Y c ~ exp[-{b x - b 2 ) 2 /r 2 ]exp[-(b 2 - hf/r% 

In the limit of r 2 ,r 2 <C R\, where Ra is a nucleus radius, the yields coincide. Thus, we 
can save the results of Ref. [E^] considering them as a summation of the yields of the 
quark-gluon diagrams. 

Let us note that neither Y a , nor Y c depend on the longitudinal coordinates or on 
the multiplicity of produced particles. It is the main difference between the "reggeon 
cascading" and the "usual" cascading. 

As well known, the intranuclear cascade model ([]T|-|§]) assumes that in a hadron- 
nucleus collision the secondary particles are produced due to an inelastic interaction of 
the projectile particle with a target nucleon. The produced particles can interact with 
other target nucleons. A distribution on distance I between the first interaction and the 
second one has a form 

11 H 

W(l) dl — exp( —I), 

</></> 

where < / >= 1/apA, cr is a hadron-nucleon cross-section, n is the multiplicity of the 
produced particles and pa — 0.15/m~ 3 is the nuclear density. At the same time the 
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amplitudes or the cross-sections of the processes shown in Fig.[| have no dependence on 
I or n. Thus, we expect that in the quark-gluon or reggeon approach the "cascade" will 
be more restricted than in the cascade model. The difference between approaches can 
lead to the different predictions for the light nuclei destruction (an effect of the limited 
volume) and for the characteristics of the heavy nuclei interactions (an influence of a large 
multiplicity of the produced particles). 

To show this, we use a simple method to estimate the nuclear destruction in the 
framework of the quark-gluon approach. 

f . As it was said above, the "reggeon cascade" is developed in the space of the impact 
parameter. Thus, for its description it is needed to determinate a probability to 
involve a nucleon into the "cascade". It is obvious that the probability depends 
on a difference of the impact coordinates of the new and the previously involved 
nucleons. Looking at the yield of the diagram of Fig. [5]c, we choose the functional 
form of the probability as 

P(\ h - % |) = C nd exp(-(b t - b,fjr n d 2 ). (15) 

Here b{ and bj are projections of the radiuses of i th and j th nucleons on the impact 
parameter plane. 

2. The "cascade" is initiated by the primary involved, wounded nucleons. If the con- 
stant C n d is small, we can use the Glauber theory for their determination. 

3. We assume that all the involved and wounded nucleons are ejected from the nucleus. 

The "cascade" looks as follows: a projectile particle interacts with some of the intranu- 
clear nucleons. They are called "wounded" nucleons. The wounded nucleons initiate the 
"cascade". A wounded nucleon can involve a spectator nucleon into the "cascade" with 
the probability fllSl). The latter can involve a second nucleon. The second nucleon can 
involve a third one, and so on. 

A Monte Carlo algorithm for estimation of the nuclear destruction in the nucleus- 
nucleus interactions which corresponds to the model formulation, includes the following 
steps: 

1. The calculation of the impact parameter distribution in the framework of the Glauber 
theory |19| ; 



2. The sampling of the impact parameter and the nucleon coordinates; 

3. The determination of the wounded nucleons (see Ref. HT9|); 



4. The determination of the spectator nucleons involved in the "cascade" by the 
wounded nucleons. If the number of the involved nucleons is equal to zero - exit; 

5. If the number of the involved nucleons is not equal to zero, a possibility is considered 
to involve the other spectators nucleons by the involved ones. If the number of the 
new involved nucleons is equal to zero - exit. In other case - it is needed to repeat 
the step 5 taking into account only the new involved nucleons. 
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The first step is performed only once at the given mass numbers of the projectile and 
target nuclei. The steps 2-5 are repeated until the needed statistics is reached. The 
steps 4, 5 are applied to the nucleons of projectile and target nuclei. 

It is suggested that all the newly involved participants and the "wounded" nucleons 
are knocked-out from the nucleus. 



3.2 Fermi-motion of nucleons 

To take into account the energy-momentum conservation law by simulating compound 
system, let us consider a reaction of the compound system (1, 2) with hadron h: (1, 2) + 
h — > 1 + 2 + h. Neglecting the transverse momenta, a final state of the reaction will be 
fully characterized by a value of merely one independent kinematical variable. As the 
variable, let us take 

xf = (E 1 + p l )/(E l + E 2 + Pl + p 2 ). 

It is useful to introduce the analogous quantity xf which satisfies obviously the relation 
xf + x~2 — 1. At given xf and xf , the rest of the kinematical variables can be determined 
by the energy-momentum conservation law. 

In case of dissociation of two compound systems A and B containing A and B con- 
stituents respectively, let us introduce for the i-th constituent of system A 

xf = (E Ai +p iz )/WX and p i± , 
and for the j-th constituent of system B 

yj = (E Bj - q jz )/Wz and q i± , 
where, E^iEsf) and Pi(q}) are energy and momentum of i-th constituent from A (B), 

W\ = J2(E Ai + p lz ), WE = J2(E Bl - q iz ). 

i=l i=l 

One can find W\ and using the energy-momentum conservation low at given 



{ x f i Pa}, {Ui i Qii}, an d determine all kinematical variables for all constituents pD |. 

According to the experimental evidence |55|, the average transverse momentum of 
spectator fragments obeys the parabolic law: 

A(A — F) 



< Pl >= -A— L < p l >; J< p l > = o.07 GeV/c. 

To reproduce this result, the values of pi± for knocked-out nucleons are simulated 
according to the distribution 



dW oc exp{-p\J < p\ >)d 2 p i± , ^<p\> = 0.07. (16) 

The sum of the transverse momenta (with sign "minus") was ascribed for the residual 
nucleus. 

The choice of xf is carried out by 

dW oc exp[-(xf - l/A) 2 /(d x /A) 2 ]dxf, < xf < 1, d x = 0.07. (17) 
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x + of the residual nucleus is determined as 1 — x t ■ 

It was assumed that the knocked-out nucleons changed their characteristics again. 
The new values of xf and pi± were simulated using the distributions (|T6| ) and ( |i7|) at 
< p 2 ± >= 0.385 (GeV/c) 2 and d x = 0.2. The results from |36], |3^, [38|] were used for 
determination of the fitting parameters. 

3.3 Excitation energy of residual nucleus. 

The change of the nuclear destruction mechanism requires a change of the nuclear residual 
excitation energy calculation procedure. In a self-consistent diagram approach one needs 
to consider the more complete diagrams and describes bound states of the nucleons. The 
last problem is not solved in the presented quark-gluon approach which allows one only to 
calculate an repulsive part of the NN-potential. It is expected that a taking into account 
the terms with high order vertex constant allows one to calculate the attractive part of 
the potential too. Until this is not made, we are to use a phenomenological approach. 
Here we follow the Ref. . 



In paper [3{| proton-nucleus interactions at intermediate energies were analyzed. The 
first stage of the interactions was considered within the Glauber approach. It was assumed 
that the projectile hadron undergoes successive interactions (elastic and inelastic) with 
target nucleon. In each of the collisions a part of the energy of the projectile hadron, E, 
is transferred to the target nucleon. The distribution on the transfer energy was chosen 
in the form: 

Fl (E) = ^e- EI{E) . (18) 

The authors of Ref. [^] supposed that the excitation energy of residual nucleus was 
the sum of the recoil nucleon energies. As a result, they described experimental data on 
neutron multiplicity dependence upon the excitation energy of the residual nuclei. 

To apply this approximation to AA-interactions at high energies, one needs to evaluate 
the number of re-scatterings of each of the knocked-out nucleons. Taking into account 
that most of AA-interactions are of peripheral nature, we assume that the slow nucleons 
knocked-out from the peripheral parts of nucleus can not penetrate deeply inside the 
nucleus because of the large NN cross-section. Their re-scatterings thus occur in the 
nearest environment. We consider the environment as the spectator nucleons being inside 
a sphere with radius r = 2 fm surrounding a wounded or involved nucleon in its initial 
state. It is assumed that each spectator nucleon may acquire recoil energy distributed 
according to Eq. (jl8|). Remember, that the nucleon coordinates were chosen randomly 
and independent according to the Saxon- Woods distribution. If a spectator nucleon was a 
neighbor of two "wounded" ones, or involved nucleons, it received the sum of two energies 
chosen according to Eq. ([Tj|), and so on. The sum of energies transferred to all spectator 
nucleons was considered as the excitation energy of the residual nucleus. 

Unlike the CEM, this method obviously will lead to the zero-excitation energy when 
all nucleons are be ejected. 

It also seems evident that the boundary between the spectator part of the nucleus 
and the part affected by the first stage of the interaction which governs the excitation 
energy of the residual nucleus, depends on the impact parameter. For example, with 
heavy projectile nuclei the excitation energy rises with decreasing impact parameter from 
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Ra + Rb to 0, approaching a maximum and then falls. According to the CEM, it must 
gradually increase. 

We have used (E) = 8 MeV and the standard evaporation model |4(J (see, also Ref. 
HI) for simulation of the residual nucleus de-excitation. 

4 Description of nucleus-nucleus interactions 

Fig. P shows the proton and 7r~-meson rapidity distributions in different nucleus-nucleus 
interactions. They were calculated by the modified FRITIOF code and by the code of 
the cascade-evaporation model^ ||41|| . Events with at least one inelastic NN-collision were 
selected at the simulations. As seen, for (id-interactions the models give close results. 
Since in this case we can neglect the cascade interactions, the coincidence of the results 
tells us about correctness of the NN-interactions description in the CEM (the CEM code 
used by us does not allow a direct simulation of NN-interactions). We can mark only a 
little enhanced baryon production in the central rapidity region in the CEM. 

In aa-interactions where the influence of the cascade interactions is sufficiently large, 
we observe a difference of the predicted 7r~-meson spectra. As expected, taking into 
account the resonances in the FRITIOF model leads to a decreased meson yield. The 
enhanced meson production in CEM becomes more pronounced in CC-interactions. 

To describe the proton production in the nuclear fragmentation regions in the frame- 
work of the FRITIOF model, we take into account both inelastic interactions of the 
nucleons considered above and elastic re-scatterings. As seen, we reproduce the baryon 
yields in the fragmentation regions at y ~ 0, 2.2 for the dd- and aa-interactions. The 
CEM calculations for CC-collisions were used for determination of the nuclear destruc- 
tion model parameters, C n d = 1, r 2 nd = 1.4 (fm 2 ). Then we introduced the experimental 
criteria for the proton registration. 

The calculations for heavy nuclei look most interesting. In Fig. [?| the proton dis- 
tributions are presented on total and transverse momentum in n + Ta interactions at 
P n = 4.2 GeV/c. As seen, the model predictions are close to each other. The same 
closeness we observe for C + Ta interactions at P = 4.2 GeV/c/nucleon presented in 
Fig. [5]. Though, there is a wide difference between the predicted meson spectra. The 
FRITIOF model calculations are close to the experimental data. At the calculation we 
use the following values of the parameters C n( i = 0.2, r 2 d = 1.1 (fm 2 ) to reproduce the 
Ta nuclei destruction. 

Summing up, we can conclude that we have reached a satisfactory description of the 
meson and nucleon production in the nucleus- nucleus interactions at the energy of 3.3 
GeV/nucleon in the framework of the sufficiently simple FRITIOF model. The model 
can be applied for practical calculation of nucleus-nucleus interaction characteristics. 

The authors of the paper are thankful to Prof. V.S. Barashenkov, Zh.Zh. Musul- 
manbekov and B.F. Kostenko for useful discussions. One of the authors (V.V.U.) thanks 
RFBR (grand No. 00-01-00307) for its financial support. 



3 'The model takes into account the trailing effect, Pauli principle, the dependence of the Fermi 
momentum on the local nuclear density, the pre-equilibrium emission and the evaporation of the nuclei. 
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Figure 1: 7r~-meson rapidity distributions in np-interactions. Histograms are the exper- 
imental data. The dashed and solid curves are the standard and modified FRITIOF 
calculations, respectively. 
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Figure 2: 7r -meson transverse momentum distributions in np-interactions. Notations are 
the same as for Fig. 1. 
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Figure 3: Proton rapidity distributions in a) np — > ppir~ , b) np — >■ ppn~7i°, c) np — >■ 
npn + 7T~ reactions at P n = 3.83 GeV/c. Histograms are the experimental data. The 
dashed and solid curves are the standard and modified FRITIOF calculations, respectively. 
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Figure 6: Proton and 7r~-meson rapidity distributions in nucleus-nucleus interactions at 
the energy of 3.3 GeV/nucleon. Points are the experimental data [42]. The dashed and 
solid curves are CEM and the modified FRITIOF calculations, respectively. 
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Figure 7: Proton rapidity distributions in nTa interactions. Points are the experimental 
data [43]. The dashed and solid curves are CEM and the modified FRITIOF calculations, 
respectively. 
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Figure 8: Proton and 7r~-meson rapidity distributions in CTa-interactions. Points are the 
experimental data [44]. The dashed and solid curves are CEM and the modified FRITIOF 
calculations, respectively. 
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